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ABSTRACT 

We have monitored the millimeter and submillimeter emission from the young 
stellar object GG Tau, a T Tauri binary system surrounded by a massive cir- 
cumbinary disk. We find that between 1992 and 1994, the flux has increased 
significantly at 800, 1100, and 1300 fim, resulting in a steepening of the observed 
spectral energy distribution at those wavelengths. Such an increase appears con- 
sistent with a modest increase in disk luminosity (a factor of two). The increase 
in the effective disk temperature might arise from a slight change in the disk 
heating processes. Alternatively, the flux increase may reflect a sudden change 
in the underlying dust optical properties. 



Subject headings: radio continuum: stars - stars: flare - stars: formation - stars: 
individual: GG Tauri 
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1. Introduction 

GG Tau is a young stellar system within the Taurus molecular cloud complex which 
is located at a distance of ~140pc (Elias 1978). It consists of a pair of close binary stars 
(Leinert et al. 1991) separated by ~10" (~1500 AU). The main component of this binary 
pair consists of a K7-M0 star (Strom et al. 1989; Hartigan et al. 1994) and an M4 star 
(Roddier et al. 1996) separated by 0.255" (~38 AU), and has collectively been classified 
as an emission- line ("classical") T Tauri star (Herbig & Bell 1988). Associated with this 
system, which we will refer to simply as GG Tau, is a large, massive circumstellar disk, 
which has recently been imaged using millimeterwave interferometers by a number of groups 
(Simon & Guilloteau 1992; Kawabe et al 1993; Koerner, Sargent & Beckwith 1993; Dutrey, 
Guilloteau & Simon 1994). They found a well- resolved, extended feature, whose morphology 
and kinematics are consistent with a rotating circumbinary disk (Dutrey, Guilloteau & Simon 
1994; Koerner, Sargent & Beckwith 1993). Dutrey et al. have also found that the disk, of 
outer radius ~800 AU, appears to have a central hole or cavity of radius ~180 AU. Roddier 
et al. (1996) have recently presented near-infrared images acquired using adaptive-optics 
techniques, which clearly show a partially evacuated cavity of radius ~200 AU in the disk. 
They also find some infrared excess from both of the stars inside the cavity, perhaps from 
warm inner disks surrounding each star. 

The circumstellar /circumbinary disk radiates very strongly at millimeter and sub- millimeter 
wavelengths (Beckwith et al 1990; Beckwith & Sargent 1991; Moriarty-Schieven et al 1994). 
Indeed, the published 800 \im flux densities are nearly 25% of the IRAS 100 /im flux density 
of 5.2 Jy. We present here evidence that the sub-mm flux density of GG Tau has varied. We 
interpret our observations as a sudden brightening or "flare" of the disk. We will present a 
set of simple models that can match the observations, and discuss their implications. 

2. Observations 

The data to be discussed in this report are presented in Table 1, and include new obser- 
vations as well as observations of Beckwith et al (1990; hereafter BSCG), Beckwith & Sargent 
(1991; hereafter BS), and Moriarty-Schieven et al (1994; hereafter MSWKT). The new ob- 
servations presented here were obtained at the James Clerk Maxwell Telescope (JCMT)Q 
located near the summit of Mauna Kea, in January, September, and December 1994. The 



1 The JCMT is operated by the Royal Observatory Edinburgh on behalf of the Particle Physics and 
Astronomy Research Council of the United Kingdom, the Netherlands Organization for Scientific Research, 
and the National Research Council of Canada. 
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Table 1. GG Tau Flux Densities 



A (/im) 


Date of Obs. 


Beam FWHM 


Su (Jy) 


References 


1 9 


1 JOO 


9^" 


i 07 _i_ n i q 


trac par* 


25 


1983 


25" 


1.65 ± 0.15 


IRAS PSC 


60 


1983 


60" 


3.03 ± 0.27 


IRAS PSC 


i nn 

1UU 




i nn" 

1UU 






450 


Jan 1994 


17.5" 


4.54 ± 0.35 


MSB 




Dec 1994 


17.5" 


4.16 ± 0.23 


MSB 


800 


29 Nov, 3 Dec 1989, 7 Dec 1990 


17" 


1.25 ± 0.08 


BS 






1 fi 'V' 
±\j.o 


i ii _i_ n i o 


MSWKT 




Jan 1994 


16.5" 


1.71 ± 0.04 


MSB 




Sept 1994 


16.5" 


1.59 ± 0.06 


MSB 




Dec 1994 


16.5" 


1.65 ± 0.08 


MSB 


i 1 nn 


9 Dpp 1 QSQ 

JJGL 1 JO J 


99" 


n onn _i_ n nci 


BS 




Sept 1991 


18.5" 


0.74 ± 0.12 


MSWKT 




Jan 1994 


18.5 


1.07 ± 0.03 


MSB 




Sept 1994 


18.5" 


0.83 ± 0.03 


MSB 




Dec 1994 


18.5" 


0.85 ± 0.03 


MSB 


1300 


29, 30 April 1988 


11" 


0.593 ± 0.053 


BSCG 




Jan 1994 


19" 


0.69 ± 0.03 


MSB 




Sept 1994 


19" 


0.63 ± 0.02 


MSB 




Dec 1994 


19" 


0.63 ± 0.03 


MSB 


2000 


Jan 1994 


28" 


0.32 ± 0.06 


MSB 



References. - 
(1991); BSCG = 
et al (1994) 



IRAS PSC = IRAS Point Source Catalog v. 2.1; BS = Beckwith & Sargent 
Beckwith et al (1990); MSB = this work; MSWKT = Moriarty-Schieven 
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photometry was obtained using the UKT14 bolometer equipped with filters passing effective 
wavelengths of 444, 790, 1090, 1260 and 1920 /im, using a 65mm aperture (corresponding to 
17.5", 16.5", 18.5", 19" and 28" (FWHM) respectively on the sky), chopping in azimuth by 
60" at a frequency of 7.8 Hz. 

No planets were available for the run of January 1994, so both atmospheric opacity 
and flux calibration were estimated from observations of CRL 618. Pointing, which varied 
by <2-3", was checked at every calibration observation performed at least once per hour. 
Flux densities were also checked using L1551-IRS5. Assumed flux densities for these two 
objects were from Sandell (1994), and are shown in Table 2Q. The quoted errors for these 
observations (Table 1) include statistical and calibration uncertainties, the latter of which 
includes the uncertainty in the flux density of CRL 618. 

During the run of September 1994, the zenith optical depth at 225 GHz as measured by 
the radiometer operated by the nearby Caltech Submillimeter Observatory (CSO), was found 
to be extremely stable, with maximum excursions ranging from r = 0.11 to 0.12. Thus we 
were able to point, focus, and obtain photometry of L1551-IRS5 (only a few degrees away) 
while it was at virtually the same airmass as GG Tau. We were also able to observe Mars ~2 
hours later when it too was at nearly the same airmass as GG Tau (1.03 compared to 1.00). 
These observations, assuming flux densities from Sandell (1994) for L1551-IRS5 and from the 
planetary flux density model of Wright (1976), were then applied directly to the photometric 
observations of GG Tau without having to estimate atmospheric opacity or instrumental 
sensitivity. (The differences in airmass between Mars and GG Tau would only introduce an 
error of ~3% even assuming Tsoo/^m ~ T) Consistent values of flux density were found from 
both calibrators, and quoted errors reflect both the statistical and calibration uncertainties. 

During the run of December 1994, we obtained the zenith optical depth at 225 GHz 
from the CSO radiometer and found the sky to be stable over the course of the observations. 
Two observations of Mars at airmasses of 1.5 and 1.1 allowed us to estimate the atmospheric 
optical depth at the observed wavelengths over the relevant airmass. We also had observa- 
tions of L1551 IRS 5 at virtually the same airmass as GG Tau. We used these two calibration 
methods (i.e. deriving r and sensitivity from Mars, and direct comparison with L1551-IRS5) 
to derive consistent values of flux density for GG Tau. 



2 In the September 1994 run we were able to check the flux densities of CRL618, and found them to be 
slightly higher than the 1991-determined values of Sandell (1994) at longer wavelengths. (These values are 
shown in Table 2.) This might be expected since Sandell noted that the contribution from free- free emission 
was slowly increasing at longer wavelengths. The effect of this would be to increase our (1994 January) 
estimated flux densities by ~17% at 1300/im and ~29% at 2mm. This would tend to enhance the flux 
increase we see between the 1992 and 1994 data sets. 
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Table 2. Flux Densities of Secondary Calibrators 



Source 


A (/an) 


Assumed" (Jy) 


Derived 6 (Jy) 


CRL 618 


450 


13 ± 2 






800 


4.3 ± 0.2 


4.3 ± 0.2 




1100 


2.6 ± 0.1 


2.8 ± 0.2 




1300 


2.3 ± 0.1 


2.7 ± 0.2 




2000 


1.7 ± 0.1 


2.2 ± 0.2 


L1551-IRS5 


450 


41 ± 4 






800 


5.7 ± 0.3 






1100 


2.4 ± 0.1 






1300 


1.7 ± 0.1 






2000 


0.7 ± 0.2 





References. — "Sandell (1994); b run of Sept 1994 
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Fig. 1. — The change of GG Tau's flux at 800, 1100, and 1300 jjm. Solid symbols are from 
this work and MSWKT, open symbols are from BS or BSCG. 

The data from BSCG were obtained at the Institut de Radio Astronomie Millimetrique 
(IRAM) 30m telescope at Pico Veleta in April 1988, and were calibrated for flux density using 
Mars and Uranus. The effective wavelength of the observations was 1300/im. BS observed 
using a bolometer on the CSO in late 1989 and 1990, calibrating for atmospheric opacity 
using HL Tau and for flux density on Mars in 1990 and on W3(OH), CRL618, IRC+10216, 
3C273, 3C84, and HL Tau in 1989 (no planets were available). The effective wavelengths 
of these observations were 1056 and 769/im with FWHM of 22" and 17" respectively. The 
observations from MSWKT were obtained at the JCMT in September 1991 at effective 
wavelengths of 790 and 1090 fim and FWHM of 16.5" and 18.5" respectively, calibrating 
flux density from Mars and atmospheric opacity from CRL618. 



3. Discussion 
3.1. The Data 

Figure [I] displays the flux density of GG Tau at 800, 1100 and 1300 fim as a function 
of the date of observation. (The 800/wn observation of BS is an average of observations 
made in late 1989 and late 1990, and so is placed at 1990.4.) The errors shown on Figure [I] 
include both statistical and calibration uncertainties, and are often dominated by the latter. 
(Note that these errors do not include the uncertainty of the absolute flux of Mars (about 
5%) based on the model of Wright (1976).) It is apparent that from September 1991 to 
January 1994 the flux density of GG Tau has increased by more than 50% at 800/im and by 



-7- 



>40% at 1100/im, but by only ~17% at 1300/im. From January to September 1994, there 
was a marginal decrease at 800/im, and a decrease to 1991 values at 1100 and 1300 \xm. 
It thus appears that the circumbinary disk of GG Tau experienced a significant increase in 
its measured luminosity, a "flare", between late 1991 and late 1993, affecting predominantly 
sub-millimeter wavelengths. 

One must use caution, however, before accepting such a claim. Errors in calibration, 
etc., could account for apparent changes in measured flux density. As noted in §2, all of these 
data except for the run of January 1994 and of December 1989 (from BS) were calibrated 
for flux density from Mars, and by using the assumed Martian flux densities calculated from 
the model of Wright (1976). The December 1989 observation at 1100/im was calibrated from 
several sources. The January 1994 run was calibrated on CRL618, whose flux density was 
checked on Mars during the September 1994 run. Thus, with minor exceptions, all of the 
calibration has been done in a consistent manner, and we are confident that the derived flux 
densities are correct within errors. 

Five of the sources including GG Tau which were observed in September 1991 by 
MSWKT were re-observed at 1100/im during our January 1994 run, providing a consis- 
tency check. These data are shown in Table 3. In four of the five cases, the ratio of 1994 
to 1991 flux densities is well within one standard deviation of being unity. Only GG Tau 
has an elevated flux density. Finally we note that the 800/im and 1100/im flux densities 
measured by BS and BSCG in 1989/1990 on different telescopes using different instruments 
are consistent with our 1991 observations and inconsistent with those of post-1991. The 
ratio of our September 1994 800/im data to the average BS 1989 data is 1.37±0.11, and the 
corresponding 1100/im data is 1.34±0.12. 

Other effects could potentially mimic a change in flux density. If the source is extended, 
a smaller aperture or beamwidth would detect less of the emission. However, interferometric 
observations, including those of Dutrey, Guilloteau & Simon (1994), have resolved the struc- 
ture of GG Tau, and have found it to be ^3" , i.e. much smaller than any of the beamwidths 
used. The component to the south associated with GG Tau(c) (the smaller of the two binary 
systems) is not visible at 2.7 mm in the Dutrey et al maps, indicating that there is little 
contribution to the millimeter emission from that source. In addition, the JCMT beam used 
by MSWKT is identical to ours, and would cover the same area. 

Finally, it must be noted that the observations of BS are made at slightly shorter 
wavelengths (769/im compared to 790/xm, and 1056/im compared to 1090//m). Since flux 
density goes as S u ~ \~ a , however, the effect of this should be to increase the flux density 
seen by BS, which is opposite to the effect seen. 
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1 1.5 2 2.5 3 3.5 

log A (pm) 

Fig. 2. — Spectral energy distribution (SED) of GG Tau. Open symbols represent 1983-1991 
data, solid symbols represent data obtained in 1994. 

We are thus confident that the increase in flux density of GG Tau seen in Figure [l] is 

real. 



3.2. Spectral Energy Distribution 

Figure |] displays the spectral energy distribution (SED) of GG Tau using the pre-1992 
data from Table 1, the 1994 January and December data, and the IRAS flux densities. 
MSWKT have shown that the slope of the SED at millimeter /sub- millimeter wavelengths 
is steeper for embedded sources compared to more evolved, optically visible, T Tauri stars 
(like GG Tau). Indeed GG Tau had the shallowest SED of all the sources in their study. 
From Figure |2] it can be seen that the slope of the SED appears to have changed as GG Tau 
has undergone the "flare". In Table 4 we show the slope a of the SED (from S v ~ A _a ) 
derived from each data set (BS includes a 624/im data point plus the data from BSCG). The 
SED appears to have steepened abruptly after the "flare" , and continues to be steeper than 
pre- "flare" . 
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Table 3. Sources Observed 1991 Sept. and 1994 Jan. at HOO^m 



Source 


S 1100 (1991 Sept) 


Shoo (1994 Jan) 


Ratio S1994/S1991 




mJy 


mJy 




04113+2758 


461 ± 53 


510 ± 40 


1.11 ± 0.15 


04169+2702 


281 ± 53 


230 ± 30 


0.82 ± 0.19 


GG Tau 


737 ± 117 


1070 ± 30 


1.45 ± 0.23 


04361+2547 


188 ± 27 


200 ± 20 


1.06 ± 0.19 


04365+2535 


438 ± 38 


390 ± 40 


0.89 ± 0.12 



Table 4. Spectral Energy Distributions 



Reference a 

BS + BSCG 1.27 + 0.15 

MSWKT 1.28 ± 0.59 

Jan 1994 1.73 ± 0.06 

Sept 1994 1.99 ± 0.10 

Dec 1994 1.82 ± 0.06 



References. — BS = Beckwith & Sargent (1991); BSCG = Beckwith et al (1990); MSWKT 
= Moriarty-Schieven et al (1994) 



-1 

1 1.5 2 2.5 3 3.5 

log A (yum) 

Fig. 3. — We display the Beckwith et al. (1990) disk model (solid line), Dutrey et al disk 
model (dashed line), and a truncated Beckwith et al. disk model (dotted line), designed to 
match the IRAS flux out to 50 jum. The truncated disk is optically thick at 50 /im and 15 
AU in radius. The nominal disk mass however is only 10~ 3 Mq. Symbols are the same as 
in Figure 2. 

3.3. Simple "Toy" Source Models 

Any model ought to account for a number of effects, including an increase of the 800 fim 
flux density of at least 50% over a timescale of <^2 years, with relatively little effect at longer 
wavelengths, resulting in a steepening of the SED during the flare and a more rapid return 
at longer wavelengths to "normal". Let us consider the existing disk models for GG Tau's 
continuum emission. In 1990, Beckwith et al presented a source model for GG Tau that 
consisted of a geometrically thin disk, with a large optically thick regime (out to 45 AU) and 
an overall disk size of 100 AU. Their model accounts for the fluxes from 12 to 1300 fim quite 
well (Fig. ^, solid line). However, Dutrey et al (1994) reported that the disk had a large 
central 180 AU radius "hole", based on 2.6 mm interferometer observations. They presented 
an alternate model, where a cold dust ring at approximately 180 AU was the dominant 
emission source of the 2.6 mm emission. Beyond the ring, which contained some 90% of 
the disk mass, the disk continued with a very flat surface density law. In addition, they 
replaced the optically thin disk with a flared disk, whose scale height was roughly 30 AU at 
180 AU radius. They argued that the temperature profile was likely to be Tocr -0 ' 5 , because 
the flared disk would intercept a larger fraction of the stellar luminosity than spatially thin 
disk models. It effectively would act more like an envelope than a disk when intercepting the 
stellar radiation field. Artymowicz and Lubow (1995) showed how such a circumbinary disk 
might arise and actually be dynamically important for the evolution of the central binary. 
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We point out that if the Dutrey et al. idea of a dense dust ring is correct, the temperature 
profile is likely to be steeper in the dust ring (and the outer disk colder) than their model 
assumes. This occurs because the dust ring, which has a large extinction, attenuates the 
radiation field of the central stars quite strongly (Wolfire and Churchwell 1994). While the 
inner temperature would likely remain at 25 K, the temperature of the outer disk would 
fall below the 20 K that is assumed in the Dutrey et al. model. The outer disk, thus, 
could be effectively shielded from the central heating source. In addition, their model does 
not produce any substantial flux below 100/xm (Fig. 0, dashed line), implying that at least 
one additional component (either a smaller circumbinary disk or circumstellar disk) must 
be present to explain the IRAS data. As an illustration of the type of disk that might be 
present, we present in Fig. |3| (dotted line) a disk model which is simply a truncated Beckwith 
et al. model. This optically thick (at 50 fim) disk has a mass of less than 10 -3 M©. With 
an outer radius of only 15 AU, it lies well inside the estimated deprojected separation of 50 
AU reported by Dutrey et al. Roddier et al. (1996) detected a near-infrared excess toward 
both stars in the cavity region, likely due to circumstellar disks. Also, the Ha emission from 
GG Tau which led to its classification as a T Tauri star is believed to be indirect evidence 
of a circumstellar disk, as opposed to the circumbinary disk seen at millimeter wavelengths. 

Given the current large uncertainties of the disk(s) physical properties, we feel that 
it is unrealistic to construct a detailed model at this time. Therefore, we will investigate 
flat disk models. For these, we adopt a single temperature power law for the disk, and 
the radii indicated by Dutrey et al., as appropriate for the submillimeter region. We adopt 
the procedures as described in Butner, Natta, & Evans (1994) to parameterize the disk 
properties. We varied the density, density profiles, temperature radial profiles, and dust 
optical properties to get a sense of how sensitive the predicted spectral energy distributions 
are to the input parameters. We found the current spectral energy distribution could be 
fit by models within a broad range of parameter space. In Fig. [|, we present three such 
models, where the "before" and "after" models are indicated by solid and dashed lines 
respectively. In each case, we limit ourselves to changing a single input parameter. In 
Fig f|(a), our model is based on a power law with a temperature profile given by T(r)oc 
r -p. p = Q 5 ^he on iy change between the models is the disk's luminosity, or effectively 
the disk temperature at the inner radius. The second model (Fig. £|(b)) has p = 0.75, 
and again, the only change between models is an increase in the disk luminosity. Both fit 
the pre-brightening phase reasonably well. By adjusting a single parameter, the total disk 
luminosity (or equivalently the temperature at the disk inner edge), we can then produce 
the post-brightening observations. The change in the disk luminosity is modest, going from 
~0.01 to ?«0.02 Lq in both cases. The corresponding disk temperatures changes are modest, 
only a few degrees with the disks remaining very cold (10-20 K). 
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An alternative to varying the temperature and large-scale disk properties is to hold these 
constant, and instead vary the dust properties. In Fig. £|(c), we present "pre-" and "post-" 
brightening models where the dust opacity law (parameterized as oc A~^) steepens from (3 = 
0.2 to (3 = 1.0. We note that our (3 = 0.2 model, while extremely high, is in agreement with 
the estimate of f3 in Beckwith and Sargent's optically thin disks (their f3 p models). 

It is possible to produce the brightening effect in other ways of course, such as varying 
the inner radius of the dust ring. The size scales involved (the nominal model disk mass is 
a substantial fraction of 1 Mq and the radius is order 200 AU), however, would argue for 
longer timescales than we see for this event. We emphasize that the fits in Fig. f| are simple 
"toy" models intended as a starting point for detailed calculations, not as a claim justifying 
a particular mechanism. 

3.4. Possible Outburst Mechanisms 

Our "toy" models illustrate that the observations do not define the parameters suf- 
ficiently to identify physical causes for the observed brightening. The effect is dramatic 
enough, however, that it is interesting to speculate. What would cause a sudden increase in 
luminosity (or temperature across the entire disk) and/or a change in the dust properties of 
a circumstellar disk? Various possibilities suggest themselves. For example, an inner disk 
might still be accreting material (Artymowicz and Lubow, 1995). If this disk were to undergo 
an outburst, it might have increased the energy available for the outer circumbinary disk 
to intercept. To date, no convincing reports of such dramatic changes have been reported 
at visual wavelengths for GG Tau (Herbst, Herbst & Grossman 1994). A standard FU-Ori- 
type outburst with its 100 times increase in luminosity certainly has not been seen. Another 
possibility is some sort of event in the circumbinary disk itself. One case that appeals to 
dynamical events is a change in the disk scale height. The change in scale height would per- 
mit the disk to intercept more radiation, thus allowing a natural way to increase the disk's 
luminosity. The distances involved — a 10 or 15 AU change in scale height is required to 
account for a two-fold increase in total disk luminosity — are large enough to be a concern. 
The required velocities are more than 20 km/sec, and it is difficult to imagine what type of 
event would cause parts of the disk to achieve such large-scale motions. Altering the grain 
properties would require a sudden density enhancement for that mechanism to be viable. A 
sudden accretion event in the disk, or interaction with GG Tau/c, the binary system which 
lies a few arcseconds to the south (cf Dutrey et al 1994), might provide the trigger. CO(2-l) 
data presented by Koerner et al. (1993) illustrates that the gas motions in the area may be 
more complex than the CO (1-0) indicates, and the near-IR images of Roddier et al. (1996) 
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show bridge-like structures between the outer ring and inner stars. If there had been an 
accretion shock as material fell on the disk, there might be additional energy - either for 
heating the grains or altering their structure enough to affect the optical properties. 

Finally, we note that one other T Tauri star, V773 Tauri (HD 283447), is known to be 
variable at mm- wavelengths (Dutrey et al. 1996; Jensen et al. (1994); BSCG; Osterloh & 
Beckwith 1995). It is also bright at cm-wavelenghts (O'Neal et al. 1990; Feigelson et al. 
1994) and is rapidly and strongly variable at these wavelengths. Its variability and large 
brightness temperature (Phillips, Lonsdale & Feigelson 1991) suggest that the emission is 
synchrotron radiation related to a strong magnetic field. It is conceivable that the sub- 
millimeter variability of GG Tau may also be related to non-thermal processes. However, 
unlike V773 Tau, GG Tau's variability seems greatest at short wavelengths, and indeed GG 
Tau has not been detected at centimeter-wavelengths despite at least two sensitive surveys 
(Bieging, Cohen & Schwartz 1984; White, Jackson & Kundu 1993). Thus it seems most likely 
that the outburst from GG Tau is thermal in nature and likely related to the circumstellar 
disk and/or circumbinary ring. 

Whatever the explanation proposed, it will require higher resolution data, particularly 
at submillimeter wavelengths, to adequately test the models. 
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Institute of Washington through a Carnegie Fellowship, as well as support from a NASA 
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Planck-Gesellschaft Arbeitsgruppe "Staub in Sternentstehungsgebieten" for their support 
while he was on leave at Jena, Germany during the writing of the final draft of this paper. 
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Fig. 4.— 3 "Toy" Models: (a) T(r) = T i (r/r i )- p , p = 0.75 where the solid line is a model 
that fits the pre-brightening phase, and the dashed line is the same model with the disk 
luminosity increased by a factor of 2; (b) T(r) = Tj(r jr^)~ v ', p = 0.50, where the solid line is 
a model that fits the pre-brightening phase, and the dashed line is the same model with the 
disk luminosity increased by a factor of 2; (c) opacity oc where (3 is changed from 0.17 
in the pre-brightening model (solid line) to 1.0 in the post-brightening model (dashed line). 
The scaling factor for the opacity law also is changed slightly. Symbols are the same as in 
Figure 2. 



